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The orientational order in vertically agitated granular rod monolayers is investigated experimen¬ 
tally and compared quantitatively with equilibrium Monte Carlo simulations and density functional 
theory. At sufficiently high number density, short rods form a tetratic state and long rods form 
a uniaxial nematic state. The length-to-width ratio at which the order changes from tetratic to 
uniaxial is around 7.3 in both experiments and simulations. This agreement illustrates the univer¬ 
sal aspects of the ordering of rod-shaped particles across equilibrium and nonequilibrium systems. 
Moreover, the assembly of granular rods into ordered states is found to be independent of the ag¬ 
itation frequency and strength, suggesting that the detailed nature of energy injection into such a 
nonequilibrium system does not play a crucial role. 

PACS numbers: 45.70.-n, 05.70.Ln, 64.70.M- 


I. INTRODUCTION 

The ordering of anisotropic particles is a universal phe¬ 
nomenon appearing widely in nature, ranging from ther¬ 
mally driven molecules or colloids [1-4] to active parti¬ 
cles such as bacteria colonies [5], actin filaments [6, 7], 
animal groups [8-10], and living liquid crystals [11]. In 
equilibrium lyotropic systems, such as hard rods inter¬ 
acting only through excluded volume interactions, the 
transition of sufficiently anisotropic particles into various 
ordered states is entropy driven. The loss in rotational 
degrees of freedom in the ordered state is compensated 
by the gain in the translational ones [3, 4, 12]. Taking a 
two-dimensional system of hard rectangles as an exam¬ 
ple, a tetratic state with four-fold rotational symmetry 
was discovered in Monte Carlo (MC) simulations [13, 14], 
and studied theoretically with density functional theory 
(DFT) [15-17]. The number density and the length-to- 
width ratio (aspect ratio) of the particles were found to 
be the key parameters determining the ordered states 
of hard rectangles with only excluded volume interac¬ 
tions [15]. Given the ubiquity of ordering transitions 
in nature, it is important to ask how well the existing 
knowledge about such transitions in equilibrium (ther¬ 
mal) systems can be extended to nonequilibrium (ather- 
mal) systems. 

Due to the dissipative interactions between particles, 
agitated granular matter has been frequently used as a 
nonequilibrium model system for phase transitions [18- 
24]. Rich and often counterintuitive dynamical behav¬ 
ior [25] has been discovered for granular rods, including 
vortex patterns [26], collective swirling motions [27], gi¬ 
ant number fluctuations [28, 29], violation of the equipar- 
tition theorem [30], and an enhanced ordering transition 
in an effective ‘thermal’ bath of spherical particles [31]. 
Reminiscent to equilibrium systems, ordering transitions 
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of vertically agitated granular rods were investigated in 
three-dimensional (3D) and quasi-two-dimensional sys¬ 
tems. In 3D, the aspect ratio of the rods was found to 
influence the ordered states of cylindrical rods [32]. In 
quasi-two-dimensional systems, a bulk isotropic-uniaxial 
nematic (I-U) transition was observed for cylindrical rods 
with large aspect ratios [33] and an effective elastic con¬ 
stant was characterized quantitatively [34]. Particularly 
in strict monolayer systems, the shape of the rods was 
found to play an important role in determining the or¬ 
dered states: Tetratic, nematic or smectic order was 
found for cylindrical rods, tapered rods or rice particles, 
respectively [35]. Moreover, tetratic order was also found 
for tubular shaped particles and the influence of the con¬ 
tainer shape was discussed in [36]. 


Despite all of these progresses, it is still unclear to 
which extent one can draw quantitative connections be¬ 
tween systems in and out of thermodynamic equilibrium. 
More specifically, a quantitative comparison between the 
state diagram of dissipative granular rods and that of the 
corresponding equilibrium system is still lacking. This 
quantitative comparison is the purpose of the present 
work. Here we investigate experimentally the orienta¬ 
tional order in monolayers of cylindrical granular rods 
driven far from thermodynamic equilibrium, and com¬ 
pare the results to MC simulations as well as DFT of the 
analogous equilibrium system. Focusing on the bulk re¬ 
gion of the system, we detect both tetratic and uniaxial 
nematic states by varying the aspect ratio of the rods. 
We demonstrate that the aspect ratio and the number 
density of rods are the key parameters determining the 
state diagram in both systems. We find a common aspect 
ratio that separates tetratic and uniaxial nematic states 
in both experiments and MC simulations. Such an agree¬ 
ment illustrates the universal aspects of the ordering of 
rod-shaped particles. 
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FIG. 1. (Color online) Sketch of the experimental set-up. The 
closed cylindrical container of height H and radius R is driven 
sinusoidally against gravity with an electromagnetic shaker. 
The rods have a length L and a diameter D. The embedded 
image shows a close view of the detected particles. 


II. METHODS 

A. Experiments 


A sketch of the experimental set-up is shown in Fig. 1 . 
Monodisperse polyvinyl chloride (PVC) rods of diameter 
D and length L, cut from welding wires of D = 3 mm 
(aspect ratio L/D < 5) or 1.5 mm {L/D > 5), are con¬ 
fined in a cylindrical container of height H and radius 
R = 10 cm. The ratio H/D = 4/3 is chosen for both 
diameters to ensure a monolayer of particles; that is, no 
rods can cross or jump over each other. The inner surface 
of the container is covered with antistatic spray (Kontakt 
Chemie, Antistatik 100 ) to minimize electrostatic forces. 
An electromagnetic shaker (Tira TV50350) is employed 
to drive the sample sinusoidally against gravity with fre¬ 
quency / = 50Hz and peak acceleration F = 47r^/^zo/^, 
where zq is the peak vibrational amplitude and g is the 
gravitational acceleration. The acceleration is monitored 
with an accelerometer (Dytran 3035B2). We capture 
high contrast images of the rods using backlight LED illu¬ 
mination and a camera (IDT MotionScope M3) mounted 
above the container. The camera is synchronized with 
the shaker so as to capture images at a fixed phase of 
each vibration cycle. The images are subjected to an 
analysis algorithm that determines the center of mass 
Pi = {xi^Hi) and the orientation Oi G [ 0 , 7 r[ of the i-rod 
with i G [1, A]. is the angle of the main rod axis with 
respect to a fixed laboratory axis, and N is the total 
number of rods in the container. The detection rate is 
100 % for D = 3 mm and 95 % for D = 1.5 mm. 

To systematically study the collective behavior of the 
rods, we vary the global area fraction 4>g = between 
^0.3 and ^0.9, and the aspect ratio L/D between 2.0 
and 13.3. For each 4>g and L/D, we vary the peak accel¬ 
eration F with a step of 1 from 2 to 20 and back. The 
waiting time between each step is fixed at ^ 1.5 minutes. 
We repeat the whole cycle at least 3 times. 


B. Monte Carlo Simulations 

Correspondingly, we model the particles as two- 
dimensional hard rectangles of length L and width D 
that interact through excluded volume interactions. N 
of such particles are placed in a box with dimensions Lx 
and Ly along the x- and ^-axes, respectively. We use 
periodic boundary conditions along both axes and study 
the equilibrium bulk configurations by means of standard 
MC simulations [37] in the canonical ensemble. That is, 
we fix the number of particles N and the system area 
A = LxLy (the temperature is irrelevant in hard mod¬ 
els). The number of particles is similar to that in the 
experiments, N ^ 10^. We use simulation boxes with 
rectangular and square shapes. No difference has been 
found between the two geometries. 

The simulation method is as follows. In order to equi¬ 
librate the system we start at very high area fractions, 
4> = ^ 0.95. We place the particles, with their main 

axes pointing in the same direction, in a rectangular lat¬ 
tice. Next we run 10^ Monte Carlo steps (MCSs). Each 
MCS is an attempt to move and rotate all the particles 
in the system. The maximum displacement Armax and 
maximum rotation A 6 >max that each particle is allowed to 
perform in a MCS is determined such that the acceptance 
probability is 0.2. Then we remove a few randomly cho¬ 
sen particles, recalculate Armax and A^max, and start a 
new simulation. The number of removed particles is such 
that the change in area fraction is A4> ^ 0.01. In order to 
rule out metastable configurations related to the prepa¬ 
ration of the initial state, we discard simulations with 
4> ^ 0.8. When the area fraction is below that limit 
we start the proper simulation. Eor each simulation we 
first run 10^ MCSs to equilibrate the system and then 
accumulate data over 10^ MCSs. Eor selected L/D we 
have also simulated the system by increasing the number 
of particles, i.e., by adding particles instead of remov¬ 
ing them. We have found no differences between both 
methods. 


C. Density functional theory 

We use an Onsager-like DET with Parsons-Lee rescal¬ 
ing. A similar DET was previously used to analyze the 
state diagram of two-dimensional rods confined in a cir¬ 
cular cavity [38]. We are interested in the behavior of 
fluid states in which the density is spatially homogeneous. 
Hence we can write, without loss of generality, the one 
body density distribution as 

= phij), ( 1 ) 

where p is the number density and h{^) is the orienta¬ 
tional distribution function. Here 7 is the angle with 
respect to the director, /^(y) is normalized such that 

d7^(7) = 1- (2) 









We split the free energy into two parts 
^[p] = ^ici[p] + ^ex[p], 
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(3) 

where Fid is the ideal gas part and Fex is the excess part 
accounting for the excluded volume interactions. The 
ideal free energy per unit of area A is given exactly by 

d 7 pft.( 7 ) ln(Apft,( 7 ) - 1), (4) 

where [3 = 1/k^T with the Boltzmann’s constant and 
T the absolute temperature. A is the (irrelevant) thermal 
volume that we set to one. The excess part is approxi¬ 
mated by 


/ d7l / d72/l{7l)/l(72)^'exc(7l2)• 

(5) 

'^exc( 7 i 2 ) is the excluded area between two rectangles 
with relative orientation 712 : 

?^exc(7i2) = {L^ + -D^)|sin7i2| +2LF(1 + |cos7i2|), (6) 

and ' 0 ex(^) is the excess free energy per particle of a ref¬ 
erence system of hard disks at the same area fraction as 
our system of hard rectangles. The diameter of the disks 
is selected such that both disks and rectangles have the 
same area. Following Bans and Colot [39] we approxi¬ 
mate T/^ex by: 

/3?/)ex(^) = (C 2 + 1) + (^2 - 1) ln(l - $), (7) 

with C 2 = 7/3 - 4 V 3 / 7 r « 0.1280. Eq. (5) recovers the 
Onsager approximation in the low density limit. 

Finally, the grand potential is given by 

^l[p\IA = F[p\IA-^p, (8) 

with // the chemical potential. We minimize Q. with re¬ 
spect to p and h( 7 ) in order to find the equilibrium states. 
We use a standard conjugated gradient method to mini¬ 
mize the functional. We use a truncated Fourier expan¬ 
sion to describe ^( 7 ). We truncate the expansion such 
that the absolute value of the last coefficient in the ex¬ 
pansion is smaller than 10 “^. 


III. RESULTS AND DISCUSSION 

This section is organized as follows: We first introduce 
the ordered states observed in experiments and MC sim¬ 
ulations in section III A. In section IIIB, we analyze the 
influence of the container walls and the driving conditions 
in the experiments. Finally in section IIIC, we quantify 
the ordering transition threshold for various aspect ratios 
and compare the state diagrams obtained experimentally, 
via MC simulations and with DFT. 



FIG. 2. (Color online) Raw experimental images (topview) 
showing typical configurations of rods with aspect ratio 
LjD — 3.3 (a) and LjD — 10.0 (b) at high global area frac¬ 
tions. The yellow (light gray) dashed circle indicates the re¬ 
gion of interest. 


A. Ordered states 


Figure 2 shows typical snapshots of the ordered states 
obtained experimentally. Short rods (a) tend to develop 
tetratic order with two alignment directions perpendicu¬ 
lar to each other. Long rods (b) form uniaxial nematic 
order with only one preferred alignment direction. In 
both cases, the container promotes either homeotropic 
(perpendicular) or planar (parallel) anchoring of the rods 
close to the boundary. To minimize the boundary effects, 
we consider only those particles located in the central re¬ 
gion of the container, as marked in Fig. 2 . A quantitative 
justification of this region of interest (ROI) will be given 
in section IIIB. Sometimes during the experiments, es¬ 
pecially at low global area fractions, we observe regions 
with very low number density of rods (almost empty re¬ 
gions). As we are interested in the bulk behavior, we dis¬ 
card those configurations in which the “empty regions” 
and the ROI overlap. 

Figure 3 shows a direct comparison of the ordered 
states obtained in both experiments and MC simulations. 
The color coded particle configurations are reconstructed 
from granular rods in the ROI (upper panels) and from 
MC simulations (middle panels) with periodic boundary 
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FIG. 3. (Color online) Typical snapshots of tetratic (left 
column, LjD — 3.3) and uniaxial nematic (right column, 
LjD — 10.0) states, (a) and (d) are reconstructed from the 
positions and orientations of the particles detected in the cen¬ 
ter region of the container, (b) and (e) are from MC simu¬ 
lations with periodic boundary conditions. The particles are 
color coded according to their orientations, (c) and (f) show 
the orientational distribution functions ^( 7 ) of the particles 
in experiments (gray bars) and simulations (solid line). 7 is 
the angle with respect to the director. 

conditions. In the tetratic state with fourfold rotational 
symmetry (left column), the orientational distribution 
function h{p{\ where 7 is the angle with respect to the 
director n, has two peaks at 7 = 0 and 7 = 7 r /2 (c). In 
contrast, in the uniaxial nematic state (right column), 
the elongated particles are oriented on average along the 
director, yielding only one peak at 7 = 0 (f). The di¬ 
rector n is calculated as the normalized eigenvector of 
the largest eigenvalue of the tensorial order parameter 
Qol^ — Here Wo^^i is the ath Cartesian 

coordinate of the unit vector Wi = (cos 6 >i, sin 6 >i), 6^^ is 
the Kronecker delta, and (...) denotes an average over 
the rods [40, 41]. To quantify the orientational order we 
measure 

qk= f d7/i(7)cos(fc7), A: = {2,4}, (9) 

Jo 

where q 2 and are the uniaxial and tetratic order pa¬ 
rameters, respectively. In an isotropic state (no orien¬ 
tational order) q 2 and q^ vanish. In a uniaxial nematic 
state ^^2 > 0 and ^4 > 0. Finally in a tetratic state q 2 = 0 
and g 4 > 0. The states in Fig. 3 are selected such that 



FIG. 4. (Golor online) Wall-rod angular correlation function 
g 4 as a function of the rescaled distance s/R to the wall for 
various L/D. A sketch with various definitions is shown in 
the inset. The data are obtained through an average over 
all F, global area fractions 4>g and experimental runs. The 
typical error for ^4 (~ 5 x 10“^) is comparable to the size of 
the symbols tor s/R < 0.5. Note the logarithmic scale on the 
y-axis. 

q 2 and q^ are comparable in both experiments and MC 
simulations. 


B. Experiments: The influence of boundary and 
driving 

Experiments [33, 35, 42, 43] and MC simulations [41] 
have shown that the container induces a preferential 
alignment of the particles close to the wall. In order 
to facilitate the investigation in the bulk, we first need 
to characterize such an influence quantitatively. 

Following the ideas in [33], we calculate the wall-rod 
angular correlation function ^ 4 ( 5 ) = (cos[4(^t7 
where 5 is the shortest distance from the rod center to 
the container wall, the angle Ot^i quantifies the tangen¬ 
tial direction of the corresponding point on the wall (see 
inset in Fig. 4), and (...) denotes an average over all the 
particles at a distance s. Either homeotropic or planar 
alignment of the particles with respect to the wall results 
in ^4 ^ 1. In Eig. 4, is presented as a function of the 
rescaled distance to the wall s/R with a binning width of 
0.03R. Eor all aspect ratios investigated, g^ decays expo¬ 
nentially with s/R. To minimize the influence of the wall, 
we consider only those particles with s/R > 0.5 to be in 
the ROT In this region, g^ is always smaller than 0.06 
and remains in a range comparable to the experimental 
uncertainties. We characterize the state of the system 
by measuring the area fraction 4> and h{'y) in circular re¬ 
gions with radius 3L inscribed in the ROT Subsequently, 
we calculate g/c(4>) from h{'y) accumulated over all the 
regions that share the same 4>. 

Eigure 5 shows the order parameters as a function of 
<F for short rods with L/D = 3.3. It indicates an area 
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FIG. 5. (Color online) Tetratic and uniaxial q 2 order pa¬ 
rameters in the ROI as a function of the area fraction $ in 
a system of rods with aspect ratio LjD — 3.3: (a) For three 
ranges of F at / = 50 Hz and (b) for four / with qk accu¬ 
mulated over all F. The threshold is obtained through 
fits to the data (straight lines) accumulated over all F for 
/ = 50 Hz (see text for details). The order parameters are 
not exactly zero in the isotropic state due to the finite size of 
the system [40]. 


fraction above which the tetratic order parameter q^ 
grows from its initial low value, while the uniaxial order 
parameter q 2 remains low. Such a combination of q 2 and 
q/^ suggests a gradual isotropic-tetratic (I-T) transition. 
As shown in (a), the behavior of q 2 and q/^ does not de¬ 
pend on the peak vibration acceleration. This is further 
confirmed through a comparison among data obtained for 
all F in the range of 2 < F <20 and also for all aspect 
ratios investigated. As shown in (b), a variation of the 
vibration frequency / from 35 Hz to 80 Hz for LjD = 3.3 
also yields the same behavior of qk{^)‘ 

Such agreements indicate that the details of how the 
rods are effectively ‘thermalized’ in our nonequilibrium 
system are not essential in determining the ordering tran¬ 
sitions, providing us the opportunity to draw connections 
to the corresponding equilibrium systems. Accordingly, 
we accumulate the data over all F at / = 50 Hz for a more 
accurate characterization of the transition threshold <Fc. 
By fitting q^^ with a constant value in the isotropic region 
and with a straight line in the ordered state, we obtain 
as the intersection point which minimizes the stan¬ 


FIG. 6. (Golor online) Tetratic and uniaxial order parameters 
as a function of the area fraction 4> for (a) L/D — 3.3 
and (b) L/D = 10.0. The insets show the corresponding 
results from MG simulations. The experimental data is an 
accumulation over all F at / = 50 Hz. Straight lines are 
linear fits to determine the threshold $c. 


dard error. Only data with sufficient statistics (i.e., error 
bar < 0.02) and q/^ < 0.3 are chosen for the fits. 

Moreover, the height of the container is found to play 
a minor role in determining the ordering transition: A 
variation oi H/D from 4/3 to 2 leads to the same be¬ 
havior of qk- Experiments with H/D = 2 for L/D = 3.3 
and L/D = 10.0 give rise to slightly lower transition 
thresholds <Fc. More specifically, we find a decrease of 
12 % for short rods and of 5 % for long rods, which is 
in both cases within the uncertainty of the fit. In addi¬ 
tion, for a specific aspect ratio of L/D = 5.0, the same 
experiments have been performed for two different rod 
diameters. The results agree with each other within the 
error bar, suggesting that the mass of the rods does not 
play a dominating role in the ordering transition. 


C. Experiments vs. simulations and DFT 

Based on the above characterizations of the boundary 
influence, we compare the ordering transitions of gran¬ 
ular rods in the ROI to the corresponding thermal sys¬ 
tem. Figure 6 shows the averaged order parameters ob¬ 
tained in both experiments and MC simulations (insets) 
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for rods with LjD = 3.3 (a) and LjD = 10.0 (b). As 
discussed above, tetratic ordering occurs in the system 
of short rods. For long rods, both order parameters start 
to grow above suggesting a gradual I-U transition. 
Qualitatively, the agreement between experiments and 
MC simulations on the behavior of both tetratic q 4 , and 
uniaxial q 2 order parameters is remarkable for both as¬ 
pect ratios. Such similarities indicate that the ordering 
of granular rods is governed by the geometric constrain of 
non-overlapping rods, which is the only interaction con¬ 
sidered in the simulations. Quantitatively, the threshold 
= 0.66 ± 0.11 obtained experimentally for rods with 
LjD = 3.3 agrees with the one 0.65 ±0.02 obtained from 
MC simulations within the error. However, the exper¬ 
imentally obtained threshold <Fc = 0.79 ± 0.04 for rods 
with L/D = 10.0 is larger than the one obtained for the 
corresponding thermal system, 0.44 ± 0.03. 

As L/D and 4> are the key parameters determining 
the state of the system, we compare the experimental 
(nonequilibrium) results with the MC (equilibrium) sim¬ 
ulations in a state diagram shown in Fig. 7. In both 
systems short rods form a tetratic state and long rods a 
uniaxial nematic state at sufficiently high area fractions. 
The aspect ratio at which the ordered state changes from 
tetratic to uniaxial nematic agrees quantitatively. It 
is found to be (L/D)t-u ~ 7.3 ± 0.7 in both experi¬ 
ments and simulations [44]. This result agrees with pre¬ 
vious simulations in which a tetratic phase was found 
for L/D = 7 and some evidence of uniaxial ordering for 
L/D = 9 [45]. The quantitative agreement of (L/T))t-u 
across systems in and out of thermodynamic equilibrium 
illustrates the universal aspects of the ordering transi¬ 
tions. 

On the other hand, the threshold 4>c for agitated rods 
differs from that in MC simulations, indicating the non- 
universal aspects of the ordering transitions. First, the 
experimentally determined 4>c exhibit a peak around 
(L/T))t-u- In contrast, MC simulations show a mono¬ 
tonic decay with L/D. Second, <Fc measured in experi¬ 
ments deviates systematically from that obtained via MC 
simulations as L/I) grows. For the largest aspect ratio in¬ 
vestigated experimentally, I//T> = 13.3, much higher area 
fraction is required for the uniaxial state to develop. This 
difference might be attributed to the following mecha¬ 
nisms. (i) The strong fluctuations in the nonequilibrium 
steady states of granular rods may lead to temporal dis¬ 
order in a system that could in principle relax into an 
ordered state, (ii) Due to the dissipative rod-rod inter¬ 
actions, the tendency of clustering for granular rods is 
larger than in MC simulations, especially for large L/D 
[compare panels (d) and (e) of Fig. 3]. (hi) Finally, the 
container wall may frustrate the orientational order of 
the agitated rods in the entire cavity. Further experi¬ 
ments using containers with different sizes and shapes 
might shed light on such a discrepancy. 

Concerning the fluctuations, it is known that the ve¬ 
locity distributions of agitated granular spheres are non- 
gaussian and exhibit exponential tails, no matter whether 



LjD 

FIG. 7. (Color online) State diagram in the plane of aspect ra¬ 
tio T/D, and area fraction 4> obtained via experiments of agi¬ 
tated granular rods (circles) and MC simulations of thermally 
driven hard rectangles (triangles). The labels denote the 
states: isotropic (I), uniaxial nematic (U), and tetratic (T). 
They are colored according to the experimental data. The di¬ 
ameter of the rods used in the experiments is D = 3.0 mm for 
L/D < 5 (gray symbols) and D = 1.5 mm for L/D > ^ (black 
symbols). Closed and open symbols indicate the I-T and I-U 
transitions, respectively. The inset shows the state diagram 
of equilibrium hard rods according to DFT in comparison to 
MC simulations in an extended region of L/D. Dashed lines 
are continuous transitions and solid lines denote first order 
transitions. 


the particles form clusters [46] or not [47]. As the dissi¬ 
pative nature does not depend on the shape of the par¬ 
ticles, we expect a similar behavior in our system. This 
feature sets agitated granular rods apart from thermally 
driven liquid crystals, and triggers the question of how to 
define an effective ‘thermal’ energy scale for a nonequi¬ 
librium system. Monitoring the mobility of individual 
granular rods with high speed photography could help 
to shed light on the difference between equilibrium and 
nonequilibrium systems found here. 

In the inset of Fig. 7 we show the state diagram accord¬ 
ing to DFT together with the thresholds obtained from 
MC simulations in an extended region of L/D. It is simi¬ 
lar to the one predicted by the scaled particle theory [15]. 
DFT also predicts I-T transitions for small L/D and I-U 
transitions for large L/D. However, the tetratic state is 
stable only for L/D ^ 2.2, most likely because only two- 
body correlations are considered in the theory [16, 45]. 
Concerning the ordering transition threshold <Fc, there is 
a good agreement between DFT and MC simulations for 
L/D ^ 7. For low aspect ratios, the deviations between 
both approaches are due to the mean field character of 
the theory. For L/D < (I//D)t-u, DFT predicts a T-U 
transition at very high area fractions. Due to the limita¬ 
tions in both experiments and MC simulations, the region 
of very high area fractions, where the T-U transition may 
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arise, has not been explored. 

IV. CONCLUSIONS 

To conclude, the ordering in agitated granular rod 
monolayers is found to be determined predominately by 
the aspect ratio of the rods and the area fraction, while 
the frequency and the strength of the agitation are not 
essential. It suggests that the detailed nature of energy 
injection into such a nonequilibrium system is not im¬ 
portant, analogous to the role that temperature plays 
in equilibrium hard rod models. In comparison to pre¬ 
vious experimental investigations on monolayer systems, 
we have focused on the bulk region of the container and 
found both tetratic and uniaxial nematic order for cylin¬ 
drical rods. This enables a direct comparison to the 
state diagram of the corresponding equilibrium system. 
We find that, depending on whether the aspect ratio is 
smaller or larger than ^ 7.3, a gradual isotropic-tetratic 
or an isotropic-uniaxial nematic transition arises as the 
area fraction grows, in both experiments and simulations. 
This agreement suggests some degree of universality for 


the ordering of rod shaped particles across systems in 
and out of thermodynamic equilibrium. Nevertheless, 
we have also found a qualitative difference between both 
systems, namely the trend of the area fraction threshold 
at the ordering transitions. 

Further investigations will focus on characterizing the 
area fraction and velocity fluctuations of the system, in 
order to find an effective ‘thermal’ energy scale for such 
a nonequilibrium system. Moreover, a comparison to 
molecular dynamics simulations [48] with tunable rod- 
rod dissipation energy could help to elucidate how fluc¬ 
tuations influence the ordering transition threshold. 
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